Residual Symmetries for Neutrino Mixing with a Large #13 and Nearly Maximal 5d 
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The residual Z|(fc) and Z^fc) symmetries induce a direct and unique phenomenological relation 
with 9 X (= #13) expressed in terms of the other two mixing angles, 8 3 (= 812) and 9 a (= #23), and the 
Dirac CP phase Sd- Z|(fc) predicts a 9 X probability distribution centered around 3° ~ 6° with an 
uncertainty of 2° to 4° while those from Z 2 (fc) are approximately a factor of two larger. Either result 
fits the T2K, MINOS and Double Chooz measurements. Alternately a prediction for the Dirac CP 
phase Sd results in a peak at ±74° (±106°) for Z|(Jfc) or ±123° (±57°) for Z 2 {k) which is consistent 
with the latest global fit. We also give a distribution for the leptonic Jarslkog invariant J v which 
can provide further tests from measurements at T2K and NO^A. 
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Introduction - The T2K Q and MINOS @ ex- 
periments indicate a relatively large reactor angle 6 X 
for neutrino mixing. At the 90% C.L., T2K gives 
0.03(0.04) < sin 2 26 x < 0.28(0.34), with zero Dirac CP 
phase, 5d, for normal (inverted) hierarchy while MINOS 
gives 0.01 (0.026) < sin 2 26 x < 0.088 (0.150) and Double 
Chooz [| with sin 2 26>i3 = 0.085 ± 0.051 at 68% C.L.. 

Many varied theoretical efforts have been made to un- 
derstand this large 6 X . Discrete groups such as S3 Q, 
A4 Q , 54 , and the binary tetrahedral group T" 
0] have been quite popular while new possibilities are 
explored in [10j. Other efforts concentrate on pertur- 
bations from some featured zeroth-order mixing such as 
democratic [H, 12 1, bimaximal Ji , 12 , T3J], tribimaximal 
@, E, 14 1, and tetra-maximal [HI patterns. More dis- 
cussions can be found in [l6| . 

In these papers symmetries or other model assign- 
ments are employed. We will show that phenomenologi- 
cal consequences of residual symmetries Z|(fc) and Z 2 (fc) 
can be readily established predicting not only 6 X to be 
large, fitting the T2K, MINOS and Double Chooz data, 
but also Sd nearly maximal in good agreement with the 
latest global fits. This provides the first strong and di- 
rect evidence for residual symmetries. 

Residual Symmetries - The symmetry that directly de- 
termines the lepton mixing pattern need not be the same 
as the full symmetry of the fundamental lagrangian. As 
the left-handed charged lepton and neutrino reside in a 
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same SU(2)l doublet, they are governed by a common 
symmetry which must be broken. Otherwise they would 
share a same diagonalization matrix (T3, [U| , leading to 
trivial leptonic mixing. It is the residual symmetry that 
determines the mixing matrices, if indeed the mixing is 
believed to be determined by symmetry. 

It is convenient to work in the diagonal basis of 
charged leptons To completely determine the mix- 
ing matrix, a product of two Z 2 symmetries is enough 
One is the well-known \i-t symmetry 



18, 2 



21] and 



the other is Z 2 [18J which can be extended to accommo- 
date a general solar angle [22| generated by, 
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There is another residual Z 2 (fc) represented by G 
G1G3, where G3 is the matrix for jjr-T symmetry [1 
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Since fi—r symmetry is just a first order approximation, 
indicated by the experimental data [lHS] and other con- 
siderations [2l( , it has to be broken. The remaining sym- 
metry would be Z|(fc) or Z 2 (fc) but not both as they are 
not independent. However, their phenomenological con- 
sequences need not be the same as we show below. Note 
that, since the diagonal mass matrix of charged leptons 
is not degenerate, G\ and G 2 only apply to the neutrino 
sector after the full symmetry is broken down to residual 
symmetries. 

Correlation Between Mixing Angles - With a single 
Z|(fc), a correlation between the three mixing angles and 
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the Dirac CP phase can be derived. In particular, 

s s — c s s x)i c a ~ s a) 



cosdo = 
cos Sd = 



(3a) 
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for Z 2 (fc) |23| and Z 2 (fc) respectively. Note that only 
physical quantities arc involved which gives the possibil- 
ity of robust physical predictions. By implementing the 
measured values of the three mixing angles, a prediction 
of Sd can be made. Or, ([3]) can be solved for 9 X , 



sin 9 X = p 



± 
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tan26» a (tan6» s ) p , (4) 



with cd = cosSd while p = ± 1 for Z 2 (fc) or Z 2 (fc). Solu- 
tions for the ± sign within the parenthesis are equivalent 
through a redefinition (9 x ,Sd) — > (—9 x ,Sd + tt) leaving 
no effect on the measured physical quantity sin 2 9 X . This 
is also true for the overall p. The difference comes from 
the exponent p leading to a (tan^) 2 w 1/2 factor be- 
tween the Z|(/c) and Z 2 (fc) predictions. 

The main feature of Q can be seen by expanding it 
to the leading order. As the reactor angle 6 X = 6 X is 
small and the atmospheric angle 9 a = 45° + S a is nearly 
maximal, ([3]) reduces to, 
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Eqs.Q-© are general and direct. To demonstrate 
this, three examples are provided. Eq.© was first ob- 
tained in a minimal seesaw model (22j with /x-r and CP 
softly broken and Z 2 (fc) retained exactly, befitting the 
situation discussed here. A special case with k = 2, 
which constrains the mixing matrix to be trimaximal, is 
studied in Q . Even an " unphysical" bimaximal solution 
24J can be covered as a marginal example. Note that 
the first two examples are obtained in mo del- dependent 
and perturbative ways while the last one comes from a 
pure symmetry analysis. 

The ratio (O of the deviation of the reactor angle 
from zero and that of the atmospheric angle from max- 
imal is given by the solar angle and the CP phase. Its 
absolute value is a minimum when Sd equals or ir, 
\S X \ > (tan 6* s ) p | . Alternately © can be solved ex- 
actly with cos Sd = ± 1 to give an absolute lower bound 
sin^ > (tan# s ) p |c a — s a |/|c a + .s a |. An upper bound can 
also be obtained but since c a ~ s a it is larger than 1. 

Numerical Predictions - A nonzero 9 X has been con- 
sistent with global fits for several years. The first hint 



appears in [25[ at only 0.9a C.L.. It persists in all subse- 
quent glob al fits [26l43d | and increases steadily to about 



3a 31 



32( as summarized in Table HI 





sin 2 6 S (0 S ) 


sin 2 a (6 a ) 


sin 2 9 X (6 X ) 


Best Fit 


0.306 (33.6°) 


0.42 (40.4°) 


0.021 (8.3°) 


la Range 


0.291-0.324 
(32.7-34.7°) 


0.39-0.50 
(38.7-45.0°) 


0.013-0.028 
(6.6-9.6°) 



TABLE I: The global fit [32|] for the neutrino mixing angles. 



The fits can be classified into two catagories depend- 
ing on the result for the atmospheric angle 9 a which is 
persistently maximal in 27, 28l 3(J 31 1 while an appar- 
ent deviation from 45° is claimed possible in 2(| 2j| 32 1 
due to a subleading effect [33l |. 
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FIG. 1: (Color Online) Predicted distributions of 6 X . 

From (U) the distribution of 9 X can be derived by using 
asymmetric Gaussian distributions P as, 

^=//^)F( S D^|, (6) 

where fg = \S(9 X — arcsins^) arc 8- functions that pick 
out the predicted value s x = RHS of for 9 X , given a 
concrete input of p, s 2 , s 2 , and Sd- There is also a mirror 
contribution for negative 9 X that is not shown - hence 
the \ pref actor. We take Sd to be evenly distributed in 
[0, 2ir), as in ©, or replace it with a specific value. The 
integration ([5]) can be simulated with scattering points 
or the delta function can be converted to one for 9 S and 
the other integrals done numerically. The results are 
shown in Fig. [T] for both Z|(fc) and Z 2 (fc). We will first 
discuss the results from Z 2 (fc). 

After averaging over Sd, the probability distribution 
peaks around 3° with an asymmetric width from 2° to 
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4°. This is in sharp contrast with the distribution given 
by the previous global fit |23[, which peaks at 0°. From 
(0 we can see that 8 X is proportional to S a . With 9 a 
significantly deviating from the maximal value, the pre- 
dicted S x must increase accordingly. In the global fit 
adopted in [23| , the central value of 9 a is about 43° and 
the maximal value is well within the la range. Hence, 
there is no apparent nonzero peak in the pre dicted dis- 
tribution of 9 X . For the latest global fit [32| of 9 a the 
central value is about 40.4°, while the maximal value is 
at the edge of the la region. This significant change in 
9 a leads to a clear nonzero prediction of 9 X . As 9 S only 
contributes as an overall factor in ([5]), its deviation will 
not change the conclusion for 9 X . For example, a differ- 
ent treatment of the reactor data leads to 3.7% difference 
(30( | in tan 9 S . The best fit values of sin 2 9 X (9 X ) vary from 
0.02 (8.1°) to 0.04 (11.5°) which are still covered by our 
predictions. This is also treated in 32J but with a much 
smaller variation, approximately 20% in sin 2 9 X . 

The measured 9 X [H-Q is not independent of the Dirac 
CP phase Sd- But this does not affect the matching be- 
tween the experimental result and the theoretical pre- 
dictions. Figure [T] shows that averaging over possible 
Sd values gives a best fit value and deviation which re- 
semble those with vanishing Sd . As the experimental fit 
of 9 X depends only slightly on Sd while the theoretical 
prediction is sensitive to it, as shown in ([5]), varying Sd 
can effectively improve the matching. For example, us- 
ing Sd = 60° moves the peak to the MINOS and Double 
Chooz central values while Sd = 70° for that of T2K. 

Since > \S a \ tan(? s , the prediction for 9 X with van- 
ishing Sd is the most conservative in the sense that it 
gives the smallest prediction for 9 X . The prediction with 
Sd uniformly distributed also peaks at 3° but has an 
extended tail to higher 9 X . This is because, while Sd 
is uniformly distributed, cos^d is not. Its distribution 
varies as (shnJo) -1 which is relatively suppressed for 
small cos Sd- Thus, the most conservative region is the 
most probable one. For example, the probability for 
\cosS D \ < 0.1(0.2,0.3,0.4) is just 6% (13%, 19%, 26%) 
corresponding to Sd = 84° (78°, 73°, 66°) respectively. 
Most of the significant region lies between Sd = 0° and 
approximately Sd — 60°. Within this region, the 9 X peak 
varies from approximately 3° to around 6° and the width 
changes from roughly 2° ~ 4° to almost 4° ^8°. This is 
the region covered by MINOS result 2.9°(4.6°) < 6 X < 
8.6°(11.4°) and 5.3° < 9 X < 10.8° for Double Chooz at 
la level while T2K has 5.0°(5.8°) < 9 X < 16.0°(17.8°) 
at 90% C.L.. 

The above discussion also applies to the case of Z 2 (fc). 
The only difference is the factor of about 2 coming from 




(a) Dirac CP Phase &u 




(b) Jarlskog Invariant J v 

FIG. 2: (Color Online) Predicted distributions of (a) the 
Dirac CP phase Sd and (b) Jarlskog Invariant J„. 



the exponent p in (jlj. As 9 X < 10° can be treated as 
small perturbation, this will induce approximate factors 
of 2 in the peak location and 1/2 in its height relative 
to the predictions from Z^/c). The result is still in good 
agreement with the data and the global fits. 

This consistency between the data and our prediction 
of a large 9 X provides the first nontrivial indication of 
the viability of residual symmetries Z 2 (fc) or Z 2 (fc). The 
correlation between the mixing angles ([3]) is independent 
of the group parameter k, and is obtained in a direct way, 
making the result quite robust. 

The change in the global fit also alters our prediction 
of the Dirac CP phase Sd [Hi- As shown in Fig. 2(a) 
the most probable value of Sd is no longer maximal. 
This is also caused by the shifted central value of 9 a . 
As S a deviates further from j , maximal S d becomes less 
probable as indicated by ©. Instead it peaks around 
±74° for Z|. Notice that a mirror solution in Q can be 
obtained through (9 x ,Sd) — > (— x ,5d + 7r) generating 
another peak around ±106°. These are in perfect con- 
sistency with the indication of —74° (—110°) for inverted 
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(normal) hierarchy [311 ]. Although no concrete number 
is provided, a nonzero CP phase also appears in (32[. For 
Z 2 (fc), the predicted cos<5_d (|3b| is larger than (j3a|) by 
a factor of 2. Consequently, the peak moves to around 
±123° (±57°). 

The distribution of the lcptonic Jarlskog invariant J„ 
is shown in Fig. 2(b) These predictions can be tested 
at T2K [3 and at NOvA (H. 



Conclusions ~ Phenomcnological consequences of the 
residual Z 2 (£;) and Z 2 (fc) symmetries arc compared with 
data and global fits. Although not independent, their 
predictions arc different. A large reactor angle 8 X peak- 
ing around 3° or 6° which is consistent with T2K, MI- 
NOS and Double Chooz can be obtained and the Dirac 
CP phase S D has peaks at ±74° (±106°) or ±123° (±57°) 
in excellent agreement with the latest global fits. This 
provides the first strong and direct support for Z|(A:) and 
Z 2 (fc) as residual symmetries of neutrino mixing. Fur- 
ther confirmation may come from the measurement of 
the leptonic Jarlskog invariant J v at T2K or NO^A. 
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